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<**] . Abstract 
• i-H ■ 

. The unexpectedly rich structure of the multisoliton solutions of the KPII 

i " i ' equation has been explored by using different approaches, running from dressing 

method to twisting transformations and to the r-function formulation. All these 
y—i • approaches proved to be useful in order to display different properties of these 

^ ' solutions and their related Jost solutions. The aim of this paper is to establish 

l/") , the explicit formulae relating all these approaches. In addition some hidden 

invariance properties of these multisoliton solutions are discussed. 

t— i ■ 1 Introduction 

The Kadomtsev-Petviashvili (KP) equation in its version called KPII 
j> ! (ut - &uu Xl + u xlXlXl ) Xl = -3u X2X2 , (1.1) 

•rH . 

where u — u(x, t), x = (xi, x%) and subscripts x\, x-i and t denote partial derivatives, is 
a (2+l)-dimensional generalization of the celebrated Korteweg-de Vries (KdV) equa- 
tion. There are two inequivalent versions of the KP equations, corresponding to 
the choice for the sign in the rhs as in (jl.ip and to the alternative choice, which is 
referred to as KPI equation. The KP equations, originally derived as a model for 
small-amplitude, long-wavelength, weakly two-dimensional waves in a weakly disper- 
sive medium [1], were known to be integrable since the beginning of the 1970s [2,3], 
and can be considered as prototypical (2+l)-dimensional integrable equations. 
The KPII equation is integrable via its association to the operator 

£(x,d x ) = -d xz +d* 1 -u(x), (1.2) 

which defines the well known equation of heat conduction, or heat equation for short. 

The spectral theory of the operator (|1.2p was developed in [4-7] in the case of a 
real potential u{x) rapidly decaying at spatial infinity, which, however, is not the most 
interesting case, since the KPII equation was just proposed in [1] in order to deal with 
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two dimensional weak transverse perturbation of the one soliton solution of the KdV. 
In fact, KPII admits a one soliton solution of the form 



u(x, t) 



^»» 2 IK 

2 2 



xi + 




2 



x 2 - 2(a 3 - b 3 )t) 



(1.3) 



where a and b are real, arbitrary constants. Multisoliton solutions have also been ob- 
tained by different methods, in [8] trough Hirota method, in [9] and [10] via dressing, 
in [11] by the Wronskian technique and in [12] by using Darboux transformations. In 
addition, in contrast with the KPI equation, also non elastic and resonant scattering of 
solitons was described in [13-17]. The problem of finding the most general iV-soliton 
solution and their interactions has recently attracted a great deal of attention. In [18] a 
sort of dressing method was applied to "superimpose" iV-solitons to a generic, smooth 
and decaying background and to obtain the corresponding Jost solutions. There, the 
two soliton case was studied in details, showing that solitons can interact inelasti- 
cally and that they can be created and annihilated. By using a finite dimensional 
version of the Sato theory for the KP hierarchy [19], in a series of papers [20-24], it 
was shown that the general AT-soliton solution can be written in term of r-functions 
and their structure was studied in details, showing that they exhibit nontrivial spatial 
interaction patterns, resonances and web structures. A survey of these results is given 
in [25] and applications to shallow water waves in [25,26]. In [27] solutions correspond- 
ing to N solitons "superimposed" to a generic smooth decaying background and the 
corresponding Jost solutions were constructed by means of twisting transformations. 

A spectral theory of the heat operator (|1.2p that also includes multisolitons has to 
be built. In [28] the inverse scattering transform for a perturbed one-soliton potential 
was derived. In [29] the initial value problem for the KPII equation with data not 
decaying along a line was linearized. However, the case of JV-solitons is still open. 
In solving the analogous problem for the nonstationary Schrodinger operator, asso- 
ciated to the KPI equation, the extended resolvent approach was introduced [30]. 
Accordingly, in order to solve the spectral problem for the heat operator, when the 
potential u(x) describes N solitons, one needs to find among different procedures used 
in deriving these potentials just that one that can be exploited in building the corre- 
sponding extended resolvent. Therefore, one needs to explore these different available 
procedures and their interrelation, which is the goal of this paper. 

The paper is organized as follows. In section 2 we introduce basic notations and 
sketch some basic topics of the standard scattering problem of operator (|1.2[) for the 
case of decaying potential. In section 3 we review the multisoliton potentials con- 
structed in [18], solving a 9-problem given by a rational transformation of generic 
spectral data. In section 4 we consider the multisoliton potentials obtained in [27] by 
means of twisting transformations and we study their connection with the potentials 
given in section 3. In section 5 we derive some alternative, more symmetric repre- 
sentations of the pure solitonic potentials, and we show that they coincide with the 
representation obtained in [20-24] by using the r-function approach. In section 6 a 
symmetric representation for the Jost solutions is also obtained, showing that they 
can be obtained by a Miwa shift [31]. Finally, we discuss the invariance properties of 
the multisoliton potentials with respect to transformations of the soliton parameters. 
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2 Background theory 

The KPII can be expressed as compatibility condition of a Lax pair, [£, T] = 0, where 
C is the heat operator defined in (| 1 . 2[) and T is given by 

T(x, d x ,d t ) =d t + Ad 3 Xl - 6ud Xl - 3u Xl - 3d-*u xa . (2.1) 

Since the heat operator is not self-dual, one has to consider at the same time its dual 
C d (x,d x ) = d X2 + d Xl — u(x) and, then, introduce the Jost solution <j>(x,k) and the 
dual Jost solution &(x, k) obeying equations 

C(x,d x )$(x,k) = 0, C d (x,d x )V(x,k) = 0, (2.2) 

where k is an arbitrary complex variable, playing the role of a spectral parameter. 

It is convenient to normalize the Jost and dual Jost solutions multiplying them by 
an exponential function as follows 

X (x,k) =e lkxi+ * 2x2 <f>(x,k), £(x,k) = e- l]iXl -* 2x2 -<&(x,k). (2.3) 

These functions satisfy the differential equations 

(-d X2 +d* l -2ikd Xl - u(x)) X (x, k) = 0, (2.4a) 
(d X2 + d 2 Xl + 2ik d Xl - u(x))£{x, k) = 0, (2.4b) 

and are chosen to obey normalization conditions at k-infinity 

lim X (z,k) = l, hm^,k) = l, (2.5) 

k — >oo k — >oo 

so that, by (|2.4p . they are related to the potential u(x) of the heat equation by means 
of the relations 

u(x) = —2i lim kd Xl x(x,k) — 2i lim kd Xl £(x, k). (2-6) 

k — >oo k — >oo 

The reality of the potential u{x) 1 that we always assume here, is equivalent to the 
conjugation properties 



X(x,k)= x (x,-k), £(z,k)=£(z,-k). (2.7) 

In the case of a potential u(x) = u (x) rapidly decaying at spatial infinity, according 
to [4-7] , the main tools in building the spectral theory of the operator (|1.2p , as in the 
one dimensional case, are the integral equations, whose solutions define the related 
normalized Jost solution xo(£,k), and £,o(x, k), i.e., 



where 



X o{x, k) = 1 + J dx' G (x - x', k)u {x')xo{x', k), (2.8a) 

^o(x,k) = l + J dx' G {x' - x,k)u {x')Zo(x' ,k), (2.8b) 

G (x,k) = -^^ f da6(a(a~2k^)x 2 )e laxl - a ^- 2W)x2 1 (2.9) 
2n J 
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is the Green's function of the bare heat operator, (— d X2 + d^ 1 )Go(x, k) = §{x). 

Thanks to (|2.8[) , the functions xo and Co have the following asymptotic behaviour 
on the x-plane 

Urn X o(a; ) k) = l J hm k) = 1. (2.10) 

x — >oo x — too 

However, if the potential u{x) does not decay at spatial infinity, as it is the case 
when line soliton solutions are considered, the integral equations (|2.8|) are ill-defined 
and one needs a more general approach. A spectral theory of the KPII equation 
that also includes solitons has been investigated using the resolvent approach. In this 
framework it was possible to develop the inverse scattering transform for a solution 
describing one soliton on a generic background [28], and to study the existence of 
the (extended) resolvent for (some) multisoliton solutions [27]. The general theory 
is, to some extent, still a work in progress. In this paper, however, we consider only 
different approaches to the construction of soliton solutions and the corresponding 
Jost solutions. 



3 Multisoliton potentials via dressing method 

In [18] we used a sort of dressing method to construct a potential u{x) describing N 
solitons superimposed to a generic background potential uq(x). More precisely, we 
considered a potential u{x) with spectral data obtained by a rational transformation 
of the spectral data of a generic background u Q (x) and we transformed the problem 
of finding this u(x) into a 9-problem on the corresponding Jost solutions, which was 
solved by a dressing procedure. The rational transformation was chosen to depend on 
N pairs of real distinct parameters aj, bj (j = 1, . . . , N) and the Jost solutions $>(x, k) 
and k) were required to obey suitable analyticity properties, normalization like 
in (|2.3|> . (|2.5|) and conjugation property (|2.7|) . 

The transformed potential u(x) is given by the following formula 

u(x) = u Q (x) - 2d 2 Xi det[C + T(x% (3.1) 

where C is a real N x N constant matrix and F{x) is an N x N matrix function of x 
with entries 

XI 

Flj{x)= J d yi ^ {y,ibi)^ {y,ia 3 ), l,j = l,...,N, (3.2) 

(bi-aj)oo 

V2=X 2 

where &o(x, k) and 1> So(x,'k) are the Jost and dual Jost solutions of equations (|2.2|) 
with potential uq(x). Notice that the matrix elements of !F(x) are given in terms of 
values of the so-called Cauchy-Baker-Akhiezer function [32] at points ctj and bi with 
j, I = 1, . . . , N. In [18] we have shown that the matrix C does not need to be either 
regular, or diagonal. We also mentioned that in order to obtain a real potential it is 
enough to consider a more general situation, with complex parameters cij and bj such 
that 

a i = S 7r a (i): bj = b„ b (j), j = 1, . . . , N, (3.3) 

where the bar denotes the complex conjugate and 7r a , 7Tb are some permutations of 
the indices (with a proper modification of the constraints on the constant matrix C). 
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However, this case is essentially more complicated to be investigated and we do not 
consider it here. 

If, a posteriori, the background uq is taken to be identically zero, then the eigen- 
functions in (|3.2p are pure exponential functions, i.e., 



where for the future convenience we introduced 

Aj(x) = ajXi + a]x 2 , Bi(x) = b l x 1 + bfx 2 , j,l = 1,2,..., N, (3.4) 
Then (|3~2)) takes the form 

f ii (x) = e- B 'WA lj e A 'W, (3.5) 

with A the Cauchy matrix 

Ay = (3.6) 

a,j - bi 

Notice that in [18] the time dependence was not specified, but this simply amounts 
to taking into account the time dependence of the original Jost solutions <&q(x, k) and 
^o(x, k) fixed by the choice of the second Lax operator (|2.ip . In the case of uq(x, t) = 
this means that we have to use 

Aj(x, t) = djXi + a 2 j x 2 - 4a^t, Bi(x, t) = b t xi + bfx 2 - 46ft, (3.7) 



instead of i 

In [18] we addressed the problem of regularity of the potentials, as well as of 
their asymptotic behaviour, only for the case N = 2, i.e., for 2-soliton potentials. In 
this particular case, we were able to formulate the conditions on the 2x2 matrix 
C that guarantee the regularity of the potential u(x), and we showed that at large 
distances in the s-plane the potential decays exponentially fast except along certain 
specific rays, where it has a solitonic one dimensional behaviour of the form (|1.3[) . The 
classification of the 2-soliton potentials obtained in [18] was then successively obtained 
and generalized to the case of ./V-soliton in [24,25]. 



4 Multisoliton potentials via twisting transforma- 
tions 

4.1 Equivalence of the TV-soliton potentials derived in [18] and 
in [27] 

In this section we consider pure soliton potentials of the heat equation obtained by 
twisting transformations. The details of the construction of the twisting operators and 
of the corresponding potentials can be found in [27]. The transformation of a generic 
(smooth, decaying at infinity) background potential uq(x) into a new potential u(x), 
which describes A^-solitons "superimposed" to the background uq(x), is parameterized 
by two sets of real parameters {ai, . . . , ajv Q } and {bi, . . . , bN b }, which we assume all 
distinct, and an N a x Nf, real matrix c. In this case we allow N a and Nf, to be not 
necessarily equal, but 

N a ,N b >l, (4.1) 
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and denote N = max{iV a , Nj,}. The pure AT-soliton potential follows directly from the 
expressions derived in [27] by taking uq(x) = 0. Precisely, we have 



u{x) = -2dl 1 \ogT 1 {x), (4.2) 

with 

n(x) = det(E Na + cF(x)) = det{E Nb + F{x)c), (4.3) 

where we introduced the N a x N a (resp. Nb x Nb) identity matrix Ejf a (resp. E^,,), 
the diagonal matrices 



e A 



<*> = diag^'W}^, e- B ^ = diag^- 8 ''*)}^, (4.4) 



see p.4[) . and the Nb x N a matrix function F(x) = || ^Fij(x)\\j_l''"' N " , where nota- 
tion (|3.5p was used, so that 

F{x) = e- B ^Ae A(x \ (4.5) 

with A a Nf, x N a constant matrix with elements given in (|3.G|) . 

Let us show now that the potentials given in (|3. 1|) in terms of N x N matrices C 
and T coincide with those obtained by the twisting transformations and expressed by 
means of matrices c and J- in (|4.3p . This is obvious when N a — Nb and matrices C 
and c are nonsingular. Then, it is enough to put c = C _1 so that the determinants 
in (|3.ip and (|4.3p are equal up to an unessential constant factor. 

The generic situation is a bit more complicated. In order to make explicit the 
size of the matrices involved we use here notation Ajy for a N X N matrix A. Let us 
consider, for defmiteness, the case N a < Nb = N. Let cm denote the N x N matrix 
constructed by adding Nb — N a zero rows to the matrix c. Taking into account that 
the product Tc in the second equality in (|4.3[) is an N x iV-matrix, it is clear that 
(T) n cn = J~c, where 

Fn = (e- B ) N A N (e A ) N (4.6) 

and where the parameters a,j with j = N a + 1, N a + 2, . . . , N in (e" 4 ) ^ and in Atv can 
be chosen arbitrarily, provided they are different from the parameters bi. It is well 
known that the Cauchy matrix An is invertible and that Ajy 1 = GnAnFn, where 
the matrix An is obtained from the matrix Ajv by renaming a,j the bj and vice- versa 
and where Gjy and Fn are two invertible diagonal matrices. Thus, for the second 
determinant in (14.311 wc have 



det(E Nb + Tc) = det(E N + F N c N ) 

= det(e- B Ae A GF) N det((c~ A ) N A N (e B ) N + G^c^ 1 ). 

The factor det(e~ B Ae A GF)iy does not contribute to the potential, due to the deriva- 
tive in (j4~2]) . Moreover, thanks to (j33]) and (|4~6]) . the NxN matrix e~ A{x ^Ae B{ - x) coin- 
cides with the N x N- matrix !F{x) in (|3.5p , once parameters aj and bj are exchanged. 
Therefore, relations ()4.3|) and (|4.2p indeed generate the same potential provided 



c N = G N CF N (4.7) 
and parameters aj are renamed bj and vice-versa. 
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Remark 4.1 In the case of a nonzero background potential uq{x) one can prove that 
the two potentials obtained by the two approaches are equivalent, up to an additional 
term logdet Tn(x). 

Remark 4.2 Here and below we omitted the time dependence as it can be easily 
switched on by using \3. 7jj instead of Jff.^p . 

4.2 Regularity conditions for the potential 

In view of the discussion on the necessary and sufficient conditions required to guar- 
antee the regularity of the multisoliton solutions of the equations in the hierarchy of 
the KPII equation, we recover here, in an equivalent formulation, sufficient conditions 
for regularity of multisoliton potentials of the heat equation and, then, of multisoliton 
solutions of KPII, already given in [21,25]. 

Let us consider again the second equality in (|4.3|) . We have 

det(E Nb +Fc) = J2 J2 (^(u ,l lT'' ,l , n )- ( 48 ) 



n=0 l<h<l 2 < - <l n <N b 



Here and in the following A ( f 1 ' ' ' ' ' \ n I denotes the minor of the matrix A obtained 



l\, I2, ■ ■ ■ , In j 

by selecting rows h,h, ■ ■ ■ ,l n and columns 31,321 ■ ■ ■ ,jn- The principal minor of the 
product Tc in (|4. 8[) can be written by the Binet-Cauchy formula as 



h, I2, ■ ■ ■ j In 



ti, «2, . . . , l n ) ^— 

l<Ji<32<---<] n <N, 




(Tc) 

\ "1) • • • 5 L n / , ^ . 

1<j 

xA( l . ul ?> " I ,.( ■'•< \ , ,4.9) 

\3x,3i,---,3nJ V h,h, ■ ■ ■ ,l n J 

Taking into account that T is a iVf, x A^ a -matrix and c a N a x ^-matrix we get that 
all minors of Tc with n > min N a , Ni, are zero. 

Recalling that the term corresponding to n = in (|4.8[) is 1, i.e., greater than 
zero, we deduce that a sufficient condition for having a regular solution is that the 
real matrix c satisfies the following characterization conditions 

a( h > l2 '--- Jn )c( K'h'-'-'iA > , (4.10) 

for any 1 < n < minjA^, Nf,} and all minors, i.e., any choice of 1 < j\ < j% < • • • < 
jn < N a and 1 < l\ < 1% < . . . l n < Nf,. In order to get a nontrivial potential under 
substitution of t\ in (|4.2[) it is necessary then to impose the condition that at least 
one of inequalities in (|4.10p is strict. 

Now, since any submatrix of a Cauchy matrix A is itself a Cauchy matrix, for an 
arbitrary minor of A we have 



Ji, 32, ■ ■ ■ , 3 V j 



^1 h,l2,...,lp\ = l<i<m<p ( j 1 1 , 
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Then, we deduce that there are orders of the parameters a's and b's for which all 
minors of the matrix A are positive. For instance, if 

ai < a 2 < ■ • • < ajv a -i < a-N a < bi < b 2 < ■ ■ ■ < fr/v^-i < bjsr b , (4-12) 
the regularity conditions (|4.10p become 

c( h ' h '---> ln )>0, (4.13) 

i.e., all minors of the matrix c must be non negative. These matrices are called totally 
nonnegative matrices [34]. Obviously one obtains the same result by using the first 
equality in (|4.3[) and the Binet-Cauchy formula enables to check directly that the two 
determinants in (|4.3J) are equal. 

If we consider the multi-time multisoliton solution of the entire hierarchy related 
to the KPII equation obtained by the Miwa shift (see [31]), as reported in (|6.22[) . we 
get again an expansion of the form (14. 9p . where the exponents are now independent 
and we conclude that the regularity condition (|4.10[) is also necessary. 



5 Equivalence with the r-function representation 

5.1 Determinant representations for the potential and Jost so- 
lutions 

In [27] we derived the transformed potential u(x) (see (|4.2I and (|4.3|) ) and the Jost 
solutions of the direct and dual problems (|2.2p related to this potential. The corre- 
sponding normalized Jost solutions are expressed as 

N " N b -Bj(x) 

X (x, k) = 1 + E e A ^ Cjl ,(E Nb + ^(x)c)- 1 J— — 

3=1 l,l'=l 1 1 

N a N h -BAx) 

3,f=l 1=1 

N a N b A . r x \ 

£(x, k) = 1 - i£ £ ^—— Cjl (E Nl + T{x)c)u}e-W 
j=i M'=i aj % 

N a N b Aj ( x ) 

sl - E Y,h^ E »« +c ^h'*- Bl{x) - ( 5 - lb ) 



— 1 f— 1 J 



j,i'=i i=i 



In analogy to (|4.3j) . all relations are given in two equivalent forms. While it is enough 
to keep only one of them, we continue to consider both of them to highlight the 
symmetry property of the whole construction with respect to the numbers N a and Nb, 
which play the role of topological charges. 

In order to obtain a representation of (|4. 3|) and (|5.ip in terms of r-functions 
we need to perform some simple algebraic operations. First, thanks to the standard 
identity for the determinant of a bordered matrix, we can rewrite (|5.ip as 

X(*,k) = ^4^, Z( X ,*) = T -^, (5.2) 
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where n is given in (|4.3[) and 



-B»(x) 



ri, x (a;, k) = det 



det 



Ti^(x, k) = dot 



det 



E Nb + T{x)c 

3=1 


6* + ik 
1 




N b 

^ &J 




1 


#jV b + T{x)c 
3=1 J 


e -S,(x) 
1 


E]y a + C T{x) 


N b 


e A,(x) 

a* + i k 


1 



ik 



--B,(x) \ 



(5.3a) 



(5.3b) 



(5.3c) 



(5.3d) 



/ 



In the above formulas, in the row and column bordering matrices EN b + J~ c and 
E^ a + cJ 7 we use the subscript * to denote an index running, respectively, from 1 to 
Nf, and from 1 to N a . Let us notice that the function t\ in (|4.3[) can also be obtained 
as a limiting value, i.e., 



T\(x) = lim ri iX (x,k) = lira Tif(a;,k). 



(5.4) 



Using (|2.3[) to go back from (|5.2[) and (|5.3[) to the Jost solutions $(a:,k) and \&(x, k), 
one can check that they have poles, respectively, at k = ibi (1 = 1,..., Nb) and k = iaj 
(j = 1, . . . , N a ) with residua 



<P b[ (x) = res $(x,k), # (x) = res *(x,k), 

k— i6/ k— ia>4 



(5.5) 



which obey the relations 



N a N b 

$ bl (x) = -i^2<$>(x,ia j )c j i, ^ aj {x) = i^cji&faibt). (5.6) 

3=1 



1=1 



Notice that these equations, together with the requirement of analyticity and the 
normalization condition (|2.10p , allow one to reconstruct the normalized Jost solutions 
in (|5.2|) and (|5.3|) and, then, via (|2.6[) the potential tt(x) in (|4.2p and 
Next, we rewrite (15.31) as 



n, x (x,k) = JJ 
\z=i 



t P--Bj(x) 



&z +ik 



det 



/ (6* + i k)[e B ( a; ) + Ae A ^c] -1, 



V 



3=1 



1 







m e ^(x) 


H 







e A ( x ) + ce 



6; + i k 
1 
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r U (z,k)= (ji e " S!(x) ) 



/ e 5 ^) + Ke Mx h 1, 



dct 



V " +ik 5 * 



- n 



a, ; + i k 



det 



N b \ 

[ e -A(x) + ce - B ^)A](a* + i k) ^ cje- 8 '^ 

1* 1 / 



By elementary transformations of the matrices on the right-hand sides, we can reduce 
them to a form where rows and columns with all elements equal to 1 or —1 are 
transformed into rows and columns with all elements equal to except 1 at the last 
place. 

In order to present the result of these transformations we introduce 



'N b 



T 2 , x (x,k) = m^ (l) (H<k) T liX (x,k), 



< x (^k) 

T2,i(x,k) 



\l=l 

No ,-Ajix) 



n 



a, ■, + i k 



n, x (x,k), 



6; +ik^ 

^,k)=(n^ (l) («i+i k ))^k). 



Then, 



JV, 



^ a,' + i k 



r 2 , x (x,k) = dct 6 w (b l+ tk)e B ^+J2^—^e A ^ Cjll 



3 = 1 



N b 



1,1' =1 



t 2 x (x, k) = det 6 



-Aj(x) N, 



33' 



T 2j? (x,k) = det 5 U 



a j + i k 



c 3i e " 



-B,(x) 



N„ 



)(h + i k) 



N a 



3,3' =1 



CjV 



k+ik (cij - bi)(a,j + ik) 

T2 >€ (x, k) = det (aj + i k)e~ A;, ' (x) + ^ 
V l=i 

Let us now introduce the limits 



Cjijbi + ik) c _ Bi(x) 



r 2 (x)= lim (ik)- JV "T 2iX (x,k) = lim (ik) iv *T 2 . ? (x,k), 

k »oo k ' oc 

r^(x)= lim( J k) Ar »T^ x (x,k)= lim (i k^r^x, k), 



(5.7a) 

(5.7b) 

(5.7c) 
(5.7d) 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 
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that, thanks to (|5.8p , have the following explicit expressions 

(\ N >> 
N a AM \ 
Su , e Bdx) +J2 € ' _°l V , (5.9a) 

r^x) = det Sjre-W + £ C f _ • (5.9b) 

By (|5.4p we have 

r 2 (x) = (JI^J ^(*) = ( Ije-^W ] n(x). (5.10) 

Both functions in (|5.9p are equivalent, in the sense that they generate the same po- 
tential 

u(x) = -2d 2 Xl \ogr 2 (x) = ~2d 2 Xl logr^x). (5.11) 
For the functions x(x,k) and £(x,k) ( see we get by flOl) and (|577|) 

Notice that functions in (|5.8p can be obtained from functions in (|5.9p by performing 
the following substitutions 

t 2 (x) -► r 2)X (ar,k) replacing e A * (oj- + i k) , e B < -» e Si (k + i k) (5.13a) 

T L{ X ) -> r 2. x ( x ' k ) re P lacin g ->• e^(ay +ik), e Si -*• e Si (6; + i k), , (5.13b) 

T2(x) — » r 2 ,f (x,k) replacing e J ^ — , e Si — * — , (5.13c) 

a.,- + i k bi + i k 

7£(a:)-»7^(x,k) replacing e^' -» 6 3 e B ' -> 6 (5.13d) 

These rules enable to significantly shorten the list of formulas given below (see also 
Remark I6.ip . 

Remark 5.1 Notice that under the transformation x —x, N a < — ► iV& and a < — > b 
the function \ is transformed into £ and viceversa. 

5.2 Symmetric representations for the potential and the com- 
parison with the r-function approach 

Here we prove that the expression (15. lip for the potential is equivalent to that one 
obtained by means of the r-function approach in the series of papers quoted in the 
introduction and surveyed in [25]. 
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We already mentioned that the double representations for the potential and the 
Jost solutions derived above highlight the symmetric role played by the parameters 
dj and 6, . To better exploit this fact we introduce N a + Nb (real) parameters 

{ki, . . .,K Na +N b } = {ai, ■ ■ ■ ,a Na ,bi . . . ,b Nb }, (5-14) 

and, by analogy with (|3.4jl . we introduce 

K n {x)=K n x l + K 2 n x 2 , n = l,...,N a + N b , (5.15) 

so that K n (x) = A n (x), n=l,...,N a , and K n (x) = B n _ Na (x), n = N a + 1, . . . , N a + 
N b . 

According to p.7p the time dependence is taken into account simply by adding a 
term — to the rhs of (|5.15[) . 

Let d denote an N a x N b real matrix with elements given in terms of the elements 
of the matrix c as 

N b N b 

d H = Hiaj-bv)- 1 ^ Yl (6,-6,0, j = l,-..,N a , l = l,...,N b . (5.16) 

l' = l l'=l,l'jil 

Let us also introduce the constant (N a + Nb) x N b - and N a x (N a + iVb)-matrices V 
and D with, respectively, the following block structures 



T>=(*\, V'=(E Na ,-d), (5.17) 
and the constant, diagonal, real (N a + Nb) x (N a + A^)-matrix 

{N a + N b "I 
IJ (Kn-«n') -1 > n=l,...,N a + N b \ . (5.18) 

Let us prove now that with one more rescaling of T2(x) and ^(x): 

!•(»)= I II (h' -bi)\r 2 {x), (5.19a) 
\i<;<;'<jV(, y 

r'(x)=l [] (-i-ai')- 1 ] (nri^-^r 1 )^^), (5.19b) 

\l<j<j'<N a J \j = l/=l / 



we get instead of 

t(x) = -A — , t'(x) = A- '- 5.20 

n ^-^ n 

l<l<l'<N b l<j<j'<N a 

where in analogy with (|4.4p we introduced the diagonal (N a + Nb) x (iV a + A^,)-matrix 

= diag{ e ^*>}^+^ (5.21) 
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and the new constant JVj x (N a + Nj,) and (N a + Nf,) x N a matrices 

K=\\K ln \\, K ln = \{ (/sn-fo/), i = l,...,JV b , n = l > ...,JV +JV 6l 
i>=ij<jki 



(5.22) 

= ll^nj Hi ^ry = ]J ( a j'- K ™)> U = 1, . . . , iV a + iV 6 , j = 1, . . . , N a . 

(5.23) 

Notice that they have a block structure since Ki t N a +v arL d ^i'j V = 1, • • • , iVt, j, j' = 
l,...,N a ) are diagonal submatrices. We also emphasize that, since N a ,Nb > 1, as 
stated in (|4.1|) . the constant matrices /C, /C', D and 2?' are not squared and, therefore, 
the determinants cannot be decomposed into products of determinants, which would 
imply u(x) = 0. 

In order to prove (|5.20[) . let us notice that 




\ n — i \ in — i Hi-Li I I 

1,1'. 

N a N b Nb N b 

= Y,e A ^d 3l , J] (a, -6,0+ E e B ""<*>«W II ~ 

3=1 l" = l,l"^l l"' = l l»=l,l"^l 

where (|5.17|) , (|5.2ip and (|5.22p were used. All terms in the last sum vanish except for 
V" = I. Thus by dH6]) we get 




IJ (6|. - &0 a det Uum^+Y, 

l<Kl'<N b I V j=l 



JV,, 



The determinant in the rhs coincides with the determinant in (|5.9ap , which proves the 
first equality in (15.20|) . The second one is reduced to (|5.9bj) by analogy, where also 
the explicit expression for matrix 7 in (|5.18j) must be taken into account. 

By considering convenient linear operations on the rows and columns of the matri- 
ces in (|5.20l) one can get more symmetric expressions of t(x) and r'(cc) which involve 
only k's. Let us consider the first equality in (|5.20p . We can subtract the last row 
of matrix IC from all the preceding rows. Then the Z-th (I > 1) row will read as 
(bi — bN b ) Yli'Li i'jti( K n — h')- Each factor bi — 6jv 6 for £ = 1, ... , Nf, can then be ex- 
tracted from the determinant and we repeat the same procedure with the last but one 
row, and up to the second one. The first row will then have all l's, and the second row 
will have entries Ki — b\ for I — 1, . . . , N a + iVj,. Then one can easily shift ni — b\ — > ni 
by using the first row. A similar transformation can be used to transform the generic 
element of the subsequent rows into K? n . Finally, instead of (|5.20|) one obtains 

r(x) = det (V e K ^ V) , t'(x) = det (V e~ K ^i V'V (5.24) 
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where by V and V we denote the "incomplete Vandermonde matrices," i.e., the Nb x 
(N a + Nb)- and (N a + N^) x iV a -matrices given as 



V= I : ; • V : : | (5.25) 

N b -1 N b -1 J I 1 JV«-1 

"1 ■ ' • K W„+« b / \ 1 ' ■ ' K N a + N b 

Remark 5.2 In the expressions i5.20\) all objects with the exception of K n (x) are 
invariant with respect to an overall shift of all parameters a 's and b's ( or, equivalently, 
all k's) by the same constant, while in \5. 25]) this invariance is not obvious. In fact, 
following the same procedure used for transforming i5.2U\) into \5.25)) one can get 
matrices V and V ' constructed from powers of K n + z instead of n n , where z is totally 
arbitrary. 

The expressions for the potential are invariant with respect to rescaling (|5.10|) 
and (|5.19p and, consequently, we have by (|4.2[) 



u(x) = -2dl \ogr(x) = -23 2 X1 logr'(x). (5.26) 

This twofold expression for the potential follows also directly, as thanks to (|5. 10|) , (|5. 14[) , (|5.15[) 
and (EH]), 

/N a +N b \ 

r(x) = (-i)^+^-i)/2 ( Jl e K n{x) \ y {Ku KNa+Nb ) T '(x), (5.27) 

where V denotes the Vandermonde determinant 

i \ ... 1 

V(Kx,...,KN a +N b ) = det j \ 

. N a +N b -1 N a +N b -1 

= Yl (k„ — K m ). (5.28) 

l<m<n<A' a +JVi, 

Thus, we have shown the equivalence of representations (|4.2|) and (|5. 1 1|) for the 
potential with those expressed as determinants of product of three and four matrices 
given in (|5.24p . These representations coincide with the representations obtained using 
the r-function approach in [21] and studied in detail in [22-25]. Notice that the special 
block form of the matrices T> and T>' in (|5.17p . in general, is not preserved, when in 
the determinants (|5.24|) we perform a renumbering of the parameters c„, for instance 
in order to have K\ < k,i < ■ ■ ■ < Kjv a +jVi,- This problem will be studied in details in 
section [ 



5.3 Explicit representation for the tau-functions 

In order to study the behaviour of the potential and Jost solutions, which we plan to 
perform in a forthcoming publication, it is convenient to derive an explicit represen- 
tations for the determinants involved, see also [20-25]. These representations involve 
the maximal minors of matrices V and T>' (see (|5.17j) ) for which we use the simplified 
notations: 

v{n u ...,n Nb )=v^ ^ ;;;< (5.29) 
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i.e., determinant of the Nt x A^-matrix, that consists of {rii, . . . ,riN b } rows of the 
matrix T> and all its columns, and 

V'(ni,...,n Na )=V ( 2 ' Na Y (5.30) 

V ni ' n 2, ■ ■ ■ , n Na J 

i.e., determinant of the N a x 7V a -matrix that consists of all rows of the matrix V' 
and {m, . . . ,riN a } columns of this matrix. Then, by using the Binet-Cauchy formula 
for the determinant of a product of matrices and notation (|5.28j) . we can rewrite 
relations (|5.24p in the form 

— , Nb 
r(x)= Yl fn^.^Yle^, (5.31a) 

l<rii<ri2<---<n]V i ,<A'a+A r b (=1 

^(*)= E fL...,n Na T[e- K ^ x \ (5.31b) 

l<n 1 <n 2 <-<n Na <N a +N b j = l 

where 

fni ,ri2,.. -,nN b = V(n ni ,. . . , K„ jVfc )D(n 1 , . . .,n Nb ), (5.32a) 

= V(n ni , . . . , K nNa )T>' '(m, . . . ,n Na )Y[^n 3 - (5.32b) 

i=i 

From (|5.27[) it follows that the coefficients / and /' in the two expansions (|5.31b[) 
and (|5.31ap are related by the equation 

jl _ ( l^N a N b +N a (N a ~l)/2 fn u n 2 ,...,n Nh (5 33) 

ni ''"'" Na V(ki, ■ ■ ■ , n Na+Nb ) ' 

where {tlx, . . . , njsr a } and {ni,n2, . . . ,npf b } are two disjoint ordered subsets of the set 
of numbers running from 1 to N a + Nj,. 

Let us mention that in our construction the two equivalent representations for the 
t- functions in (|5.27p were obtained as a consequence of the two equalities in (|2.6[) , 
while, of course, they can be proved directly by means of (|5.31|) and (|5.33[) . Sec [24], 
where this property is called duality. 

Remark 5.3 Notice that, in agreement with and \5.33\) . we have 

N, 



3=1 

and 



IJ 7n, V( Kni , . . . , . nNa ) = det ^(-1)^+ 1^1 — ^ ) , (5.34) 

V{Kl, ■ ■ ■ ,K Na+Nb ) 



T>'(ni, . . .,njyj = det7rX»(ni, . . -,n Nb ), (5.35) 
where tt is the matrix performing the permutation from (n ni , . . . , K n]Vn , Kn t , ■ • ■ , ) 

to («i, . . . , K Na+Nb ). 
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6 Jost solutions and invariance properties 

6.1 Properties of matrices V and V' 

Matrices T> and V' introduced in (|5. 1T[> obey rather interesting properties. As follows 
directly from the definition, they are orthogonal in the sense that 

D'D = 0, (6.1) 

where zero in the rhs is a N a x ^-matrix. Moreover, since the matrices 

V^V = E Nb +d)d 1 V'V' ] = E Na +dd\ (6.2) 

where | denotes Hermitian conjugation of matrices (in fact, transposition here), are 
invertible, the matrices (see [33]) 

{V) { - 1] = {E Nb +Sd)- 1 V\ (6.3) 
=V'\E Na +dd?)-\ (6.4) 

are, respectively, the left inverse of the matrix T> and the right inverse of the matrix 
V, i.e., 

(£>) M) V = E Nb , V'(V') {1) =E Na . (6.5) 

Products of these matrices in the opposite order give the real self-adjoint (N a + Nb) x 
(N a + iVb)-matrices 

P = V(V) { ~ 1) =V(E Nb +dUy 1 (V)\ (6.6) 
P' = (V'Y^V = (V')\E Na +dd))- 1 V 1 (6.7) 

which are orthogonal projectors, i.e., 

P 2 = P, (P') 2 =P', PP' = = P'P, (6.8) 

and complementary in the sense that 

P + P' = E Na+Nb . (6.9) 

Orthogonality of the projectors follows from (|6.ip and the last equality from obvious 
relations of the kind (E Nb + d t d)~ 1 d^ = d^{E Na + dd^y 1 . 



6.2 Symmetric representations for the Jost solutions 

In order to get a r-representation for the Jost solutions, we use (|5.12p and notice that 

rescaling (|5.19|) does not modify the substitution rules given in (|5.13|) . In terms of 
notations (|5.14p . (|5.15p and (|5.2ip these rules read as 

t(x) — » t x (x, k) replacing e K ^e K (n + ik), (6.10a) 

t'(x) -»T^(x,k) replacing e K e K (n + ik), (6.10b) 

t(x) — » t^(x, k) replacing e — > — , (6.10c) 
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t'(x) — » r£(x,k) replacing e K — > — - — — . (6.10d) 

where n + ik denotes the diagonal (N a + Nb) x (N a + iV{,)-matrix 

k + i k = diag{^ 1 + i k, . . . , K Na+Nb + i k} (6-11) 
and analogously for the matrix (re+i k) -1 . Explicitly these replacements give (see (|5.24[) ) 

r x (x,k) = det (ve K(x) (K + ik)p) , (6.12a) 

T^(x,k) = det (D'e-^fc + ikJ-Vj , (6.12b) 

T{(s,k) = det (Ve^^re + ik)" 1 !)) , (6.12c) 

7^(s,k) = det fv' e- K ^(K + ik) 1 V'^j , (6.12d) 

Thus, by (j5~T2|) we have 

[ll(6i+ik)J X (*,k) = ^ s ^ n ^+* k )J^p ( 6 - i3a ) 
[fifo**)- 1 ) ^ = ^ 55 (If K +zkrl ) ^%r- (fU3b) 

The Jost solutions themselves are then given by (|2.3[) and in order to simplify their 
analyticity properties it is convenient to renormalize them in the following way 

*(x, k) - N Hx ' k) , k) - ( JJ(6, + i k)\ k). (6.14) 

II(6i + <k) W 7 

Then, to preserve relations given in (12. 3p we also normalize xC^k) and £(x,k) ac- 
cordingly, so that instead of (|6.13|) we have 

Now xO^j k) is a polynomial with respect to k of the order k^ 6 and £(x, k) is a mero- 
morphic function of k that becomes a polynomial of the order k^" after multiplication 
by n^i"^ ( K n+i k). In other words, now &(x, k) is an entire function of k and ^(x, k) 
is meromorphic with poles at all points k = k,„ n = 1, . . . , N a + iVj,. Introducing the 
discrete values of $(x,k) at these points as a (iV + iVf,)-row 

$(x, ire) = {$(x, ire x ), . . . , $(x, zre^+^J}, (6.16) 



A' 
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and the residuals of ^(a;,k) at these points 



V Kn (x) = res *(af,k), (6.17) 



as a (N a + A^)-column 



9 H (x) = Kl (x), . . . , * KNa+Nb (x)} T , (6.18) 

we get thanks to (|5.14[) . (15. 16[) and (|5 . 1T[) that the relations (|5.6p take the more 
symmetric form 

<t>(x,iK)V = 0, V'V K (x) =0. (6.19) 

It is necessary to mention that after the renormalization (|6.13[) the asymptotic condi- 
tions l|2.5p become 

lim (i k)- Nb X (x, k) = 1, lim (i k) N »£(x, k) = 1, (6.20) 

k — >oc k — >oo 

and the relations (|2.6p take the form 

u(x) = -2 lim (i\s)- Nb+1 d Xl x(x,k) = 2 lim (i k) Nb+1 d Xl £(x, k). (6.21) 

k— »oo k • oc 

Finally let us point out that by introducing an infinite set of times, and, precisely, 
replacing K n (x) with the formal series (see [31]) 

oo 

K n {t 1 ,t 2 ,...) = J2< t j, (6-22) 

appropriate choices of the times tj's provide the multisoliton solutions of any non- 
linear evolution equation in the hierarchy related to the KPII equation, as well as 
the corresponding Jost solutions. In particular the multisoliton solutions of the KPII 
equation are obtained by choosing t\ = Xi, t 2 = x 2 , t$ = —At and all highest times 
equal to zero. 

Remark 6.1 We showed in \6.10)) that the Jost solutions can be obtained by means 
of transformations, that are equivalent to the formal Miwa shift used in the construc- 
tion of the Baker- Akhiezer functions in terms of the r-functions. In fact, substitu- 
tions \6. 1 0b\) and h6.10a\) can be obtained (up to an unessential factor by con- 
sidering the multi-time r-function obtained by replacing K n with the infinite formal 
series in 16'. 22\) and, then, by shifting tj — » tj — l/j(i/k) J . Similarly, the substitu- 
tions H 6. 10a]) and \6.10c\) are obtained (again up to an unessential factor —i/\s.) by 
the shifts tj — > tj + 1 / j (i / k) J . Then, non formal Jost solutions of the heat equation 
with a potential being a solution of KPII are derived by choosing t\ — X\, t 2 — x 2 , 
t^ = —At and all highest times equal to zero. 

6.3 Invariance properties of the multisoliton potential 

In some cases it is useful to rename the spectral parameters n n ->K n , for instance in 
such a way that the renamed parameters Kji are ordered as follows 

Kl < K 2 < ■ ■ ■ < (6.23) 
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This permutation can be performed by means of a (N a + Nb) X (N a + A^-matrix tt 
such that 

(k 1} . . . , K Na+Nb ) = (/%,..., KjVa+JVb) 71 "- (6-24) 

This matrix is unitary, 

(6.25) 



7T^ = 7T 1 , 



and all rows and columns have one element equal to 1 and all other equal to 0. It is 
convenient to write this matrix in a block form like 



I'll ""12 
7T21 ""22 



(6.26) 



where n%x is a N a x A^-matrix, 7r 2 2 a A^ X A^-matrix, 7r 12 a A^ x AVmatrix and 7T2i 
a A^ x A^-matrix. 

Then, representations (|5.24[) under the transformation (|6.24[) keep the same form 
if the matrices V and V ' are transformed as follows 

D->2? = 7rZ> =(*;), V ->v' = V'ir1 = (d[,-4), (6.27) 



where, from (|5. 1 7|) . we have 

d\ = Wild + 7Tl2, d 2 — TT 2 ld + 7T 2 2, (6.28) 

dj[ = ti^ - d7r{ 2) -d 2 ' = 7r| x - d7r| 2 . (6-29) 

While the size of blocks in (I6.27|) are the same as in (|5.17|) . the block structure in 
general is different. Nevertheless, thanks to the unitarity of the matrix tt, relation (|6.ip 
remains valid for the transformed matrices, which by (|6 . 2T[) means that 

d^di = d' 2 d 2 . (6.30) 

Also relations (|6.19p are preserved for the transformed quantities and transformed 
projectors P and P 1 can be built. 

Let us mention that the special block structure of the matrices T> and T> ' in (|5.1T[) 
determines them uniquely in correspondence to a given potential. One can release 
this condition, without changing the potential, by multiplying the matrix T> by any 
nonsingular Nj, x At,-matrix from the right and the matrix T>' by any nonsingular 
N a x A^-matrix from the left. In fact, determinants of these matrices cancel out 
in (|5.26[) and (|6.13[) . Viceversa one can use this procedure for bringing T> and V 
back from the block structure (|6.27p to the special block structure in (|5.17[) . This can 
be done by using matrices (d2) _1 and (d{) _1 , if they exist, to perform the following 
transformations 



di(d 



V^V(d 2 )-' = j , V " - {d[)- L V = (E Na ,-(d[)-'d>), (6.31) 

and, then, by noticing that thanks to (|6.30[) 

d^y 1 = (df) -1 ^ =d. (6.32) 

Taking into account that both representations in (|5.24p are equivalent, we deduce 
that matrices d 2 and d[, if invertible, are simultaneously invertible. Thus, we proved 
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that in the case of a permutation it such that the matrix c?2 (or d[) is nonsingular, 
the permutation of k's is equivalent to a transformation of the matrix d to d, or, 
correspondingly, by (|5.16[) . to a transformation of the matrix c. This is always the 
case when the matrix ir has a diagonal block structure, i.e., when -k\2 = ■K21 = 0. 
In this situation both matrices di and d[ are invertible and one can use the above 
substitution. Notice that in this case the permutation of k's does not mix the original 
parameters a's with 6's (see (|5.14|Q . 

In general in making a permutation of k's, say, reducing them to the order given 
in (|6.23|) , the block structure (|5 . 1 7[) is lost and we can only say that the r-functions are 
given by (|5.24|) and (|6.12[) , where both matrices T> and V ' have at least two nonzero 
maximal minors. 

7 Concluding remarks 

Here we described relations between different representations, existing in the literature, 
of the multisoliton potentials of the heat operator and we derived forms of these 
representations that will enable, in a forthcoming publication, a detailed study of 
the asymptotic behaviour of the potentials themselves and the corresponding Jost 
solutions on the x-plane. We also presented various formulations of the conditions 
that guarantee the regularity of the potential. Nevertheless, the essential problems 
of determining the necessary conditions of regularity of the multisoliton potential is 
left open. In this context, let us mention the special interest of the specific subclass 
of potentials satisfying strict inequalities in (|4.10p , which by (j 5 . 3 1 ap is equivalent to 
requiring that all f ni , . . . ,N b are of the same sign (analogously, by (|5.31b[) that all 
/ni > • • • >iV are of the same signs). These conditions identify fully resonant soliton 
solutions (cf. [20, 21]). When such conditions are imposed, all maximal minors of 
matrices V and T> are different from zero, as follows from (|5.32a[) and (|5.32b[) and, 
then, thanks to the invariance properties discussed above, we can always permute 
parameters k's in any way, for instance as in (|6.23p . and at the same time deal with 

transformed matrices T> and T> , which have a special block structure as in (|5.17[) . 
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